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cooling rate in the strong drive regime. Further we discuss the trade-off between 
efficiency of cavity cooling and robustness with respect to ubiquitous imperfections in 
a realistic experimental setup, such as anharmonicity of the trapping potential. 
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1. Introduction 

Cold molecules are nowadays very actively explored pQ. Of particular interest are cold 
heteronuclear molecules prepared in their electronic and vibrational ground state, which 
in view of their electric dipole moments can provide a strong dipole-dipole interaction, 
and a strong coupling of rotational excitations to microwave fields and cavities. The 
field of potential applications of cold samples of polar molecules ranges over ultracold 
chemistry [2] and precision measurements 0, H] , quantum computation [5] [6] and the 
simulation of exotic condensed matter models j7]|8l|9]. A necessary prerequisite for these 
developments is the ability to cool molecular ensembles and single trapped molecules to 
very low temperatures. 

Several techniques for cooling, and preparation of ultracold molecular ensembles 
have been investigated [H)] both experimentally and theoretically, including cavity- 
assisted laser cooling [11], [121 HH], buffer gas cooling, and Stark deceleration of polar 
molecules including trapping in magnetic and electric traps [Ml [15] [16], [TT], ITS] , 
evaporative [19] and sympathetic cooling schemes [20], and preparation of ground state 
molecules by photoassociation of cold atoms [2TL 1221 1231 [24]. Recently a technique 
for trapping and cooling of polar molecules on a chip has been suggested, which 
relies on strong coupling of rotational states of the molecule to a superconducting 
stripline cavity [25] . When trapped at a distance d ~ 1 fim above the cavity 
electrodes, the vacuum Rabi frequency g associated with the exchange of a rotational 
excitation and a single microwave photon is in the order of 100 kHz. In combination 
with a comparable photon loss rate 2k ~ g this coupling serves as an efficient 
dissipation channel for (the in free space essentially stable) rotational states, which 
is exploited for cooling the motion of a trapped molecule. Cooling forces in this 
setup arise from the strong gradients of microwave fields in the near-field regime, in 
analogy to cavity-assisted laser cooling techniques discussed for atoms in the optical 
regime [11] [26] [27] [28] [29] [HQ] [HI] ■ Microwave cooling does not rely on specific electronic 
or collisional properties [2] [19] [20] [32] [33] [34] of the molecules and therefore should be 
applicable for a large set of di- and also multiatomic polar molecules. As molecules 
are prepared in the motional ground state of a chip-based trap, this cooling technique 
is in particular interesting for applications in the context of on-chip (hybrid) quantum 
information processing [25] [35j. 

In this work we will present a detailed theoretical study of cavity assisted cooling 
of a single trapped polar molecule in the microwave domain. The goal of our theoretical 
analysis is to identify the physical mechanisms behind different cooling schemes, and to 
discuss the cooling rates and final temperatures in various parameter regimes including 
weak and strong driving fields and a finite temperature of the environment. The 
main conclusions of this analysis are as follows. Depending on the detuning between 
the relevant rotational transition frequency and the microwave fields, cavity cooling 
is dominated either by cavity assisted sideband cooling (CASC), where cooling forces 
arise from the gradient of a classical driving field, or "Vg" cooling where forces are 
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due to gradients in the coupling between the molecule and a single microwave photon. 
Zero environmental temperature in principle allows ground state cooling in the resolved 
sideband limit for the respective processes. At finite temperature T, thermal cavity 
photons at frequency uo c introduce an additional heating mechanism, and the minimal 
temperature Tf is given by Tf = (v/oj c )T <C T with v the trap frequency. We find 
that this limit can be reached in the strong driving regime, where also cooling rates 
are optimized. Our analysis of cavity cooling in the presence of imperfections such as 
anharmonicity of the trapping potential shows that by switching between the CASC 
and the Vg cooling regime, cooling efficiency is traded for an increased robustness of 
the cooling scheme. This may in particular be important for ongoing experiments. 

The paper is structured as follows. In section [2] we introduce the model for a 
trapped polar molecule coupled to a stripline cavity, and present a brief overview of the 
basic physical ideas and mechanisms underlying the cavity assisted microwave cooling. 
In section [3] we proceed with a technical derivation of the cooling and heating rates in 
the bad cavity limit. Based on these results we present in section H] a more detailed 
discussion of cooling rates and final temperatures as a function of different system 
parameters. Finally, in section [5] we summarize the main results and conclusions of 
this work. 

2. Overview on cavity-assisted microwave cooling of a trapped polar 
molecule 

The goal of this work is to present a theory for ground-state cooling of trapped polar 
molecules, in analogy to the ideas and theoretical tools employed for cooling of trapped 
ions in the optical domain [36l [37] . In this section we start with a brief review of 
the system involving a polar molecule coupled to a stripline cavity, as proposed in 
[25] . Further, the general relation between cooling/heating rates and the spectrum of 
fluctuations of the dipole forces which act on the molecule allows us to identify different 
cooling mechanisms in this setup. A qualitative discussion of the individual mechanisms 
is presented together with analytical results for a weakly driven molecule. The goal of 
this overview is to provide a physical understanding of the main cooling and heating 
processes which serves as a guideline through the more technical derivation of the force 
spectrum in section [3] and the extended discussion in section HI 

2.1. Cavity QED with a single trapped polar molecule coupled to a microwave cavity 

We consider the setup shown in figure where a single polar molecule is trapped close 
to a superconducting stripline cavity. The system Hamiltonian H$, 

Hs = H c + H m + H int , (1) 

is given in terms of the Hamiltonian H c for the microwave cavity, the Hamiltonian H m 
describing the internal degrees of freedom of the molecule including a classical microwave 
driving field, as well as its motion (external degree of freedom) in the trap. Finally 
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Figure 1. a) Schematic plot of the setup considered in this paper for cavity-assisted 
cooling of polar molecules. A more realistic trap design can be found elsewhere [25] , 
b) Level diagram for the cavity-molecule system restricted to the lowest states |<?)|0) c , 
|e)|0) c and |<?)|l) c . Dashed lines indicate different quantized motional levels. See text 
for more details. 

H int describes the dipole interaction between rotational states of the molecule and the 
electric field of the cavity. Apart from the coherent evolution governed by H§, coupling 
of microwave photons to the electromagnetic environment introduces a dissipative decay 
channel for cavity photons. In section [3] we incorporate dissipation in terms of a master 
equation description of our system. In the following we summarize the properties of Hs 
which are relevant for the discussion in this paper while more details can be found in 
[251 EH]. 

With initial temperatures of the molecule in the order of 1-100 mK, electronic 
and vibrational degrees of freedom are frozen out and the molecule is well described 
by its rotational degrees of freedom only. The anharmonicity of the rotor spectrum, 
E TOt (N) = BN(N + 1), with B the rotational constant and N the angular momentum 
quantum number, allows us to identify two specific rotational states \g) and |e) with 
a transition frequency uo eg = (E e — E g )/h typically in the tens of GHz regime. With 
spontaneous emission rates of low lying rotational states in the mHz regime for zero 
temperature, the idea is to create an effective decay channel for the excited state |e) 
via a microwave stripline cavity. When the frequency of the relevant cavity mode u c is 
near resonance with u eg , the dynamics of the molecule is restricted to the two states 
{\g), |e)} and the molecule is well approximated by a two-level system (TLS). 

The microwave cavity sketched in figure [I] is a (quasi-) one-dimensional 
superconducting stripline with length L ~ 1 cm, and transverse dimension of order 
d ~ 1 fim. Allowed wavevectors k n = nn/L correspond to a mode spacing in the GHz 
regime. Restricted to a single resonator mode, the cavity Hamiltonian is H c = fiw c c^c 
with c (c^) the annihilation (creation) operators for photons of frequency u c . The electric 
field generated by the cavity is given by the operator E c = So(f)(c + c') with £o(r) the 
field strength associated with a single microwave photon. While £q{t) varies on the scale 
of L along the cavity axis, the spacial dependence of So(f) in the transverse direction 
or at the cavity end-caps is determined by the typical electrode spacing d. For strong 
cooling forces we are in particular interested in those strong gradients in the near field 
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regime. The microwave photons dissipate through the coupling of the stripline to its 
surrounding leads. At frequencies of u> c ~ 10 GHz, quality factors in the range of 
Q ~ 10 3 — 10 6 [391 SOI SI] translate into a photon loss rate Ik in the range of MHz to 
kHz. 

Molecules interact with microwave photons via the dipole interaction H mt = —fi-E c , 
where (i is the electric dipole operator of the molecule. The rotational states \g) and |e) 
are coupled to the quantized electric field of the cavity. Disregarding, for the moment, 
the motional degrees of freedom we can describe the combined cavity-molecule system, 
including a classical microwave field of frequency u m and Rabi-frequency Q which drives 
transitions between \g) and |e), by a Jaynes-Cummings type model in a frame rotating 
with the microwave frequency u m , 

HA HQ 
H JC = h (Ac - A) c f c -—a z + hg (a+c + <r„c f ) + — a x , (2) 

where A c = u c — u> eg and A = u m — u eg are the detunings of the cavity field and the 
classical driving field from the rotational transition frequency, respectively, and we have 
adopted the usual Pauli operator notation for the TLS. Under resonance conditions 
uj c pd u eg equation (j2J) describes coherent oscillations between an excitation in the TLS, 
|e)|n) c , and an additional photon in the cavity, \g)\n + l) c , on a timescale (gy/n + 
[121 US]- For a molecule trapped at a distance h ~ d ~ 1 — 0.1 fim above the cavity and a 
molecular dipole moment of the order of several Debye the vacuum Rabi frequency g is 
of the order of In x 50 — 500 kHz [25]. In combination with this coherent energy transfer 
between molecule and cavity, the photon decay opens an effective decay channel for the 
rotationally excited state \e). For k > g the corresponding characteristic decay rate is 
given by 

7= ^TAf (3) 

Although the validity of this expression is limited to 7 < g, we find a significant increase 
of the rotational decay rate from T eg ~ mHz in free space to a cavity enhanced decay 
of 7 ~ 10 — 100 kHz. This also implies that in principle cooling rates of the same order 
can be achieved in this system. 

Trapping of a molecule close to the surface of a chip requires a sufficiently strong 
trapping potential; electrostatic traps with frequencies of the order of v/2ix ~ 1 MHz 
are experimentally feasible [25] 0. While in general the trapping potentials will be 
different for the two states \g) and |e), there exist so-called 'sweet spots' for the electric 
bias fields, where this difference vanishes [25]. Let us consider a ID model in which 
case the external degrees of freedom are described by a harmonic oscillator Hamiltonian 
He = hua^a, with a (a') the usual annihilation (creation) operators. Coupling between 
internal and external degrees of freedom arises from the spatial dependence of the Rabi 
frequencies g — >• g(x) and Q — > Q(x), as compared to the case of laser cooling which we 

I In principle our derivation applies equally well to magnetic traps, as used in current atom chip 
experiments |44j , in the case of a molecule with an unpaired electron spin. Note however that the trap 
fields must be compatible with the superconducting setup. 
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already remarked in the introduction. Here x is the position operator of the molecule 
with respect to the center of the trap. Including motional degrees of freedom into 
equation (j2J), we end up with the explicit form for the system Hamiltonian Hs ([I]), 



The bare energy levels of Hg and the different couplings between them are sketched 
in figure [I] b). We emphasize at this point that the validity of if m given in equation 
(J1J) is based on two crucial assumptions: i) the trapping potential is harmonic and ii) 
the trapping potential is state independent. Assumption i) is in general fulfilled near 
the trap center, i.e., in the final stage of the cooling process and as mentioned above, 
assumption ii) is fulfilled for certain conditions on the trapping fields. As neither of 
those assumptions will be strictly satisfied in a realistic experiment, we will investigate 
small deviations from these conditions in section HI Note that as long as condition i) 
is approximately fulfilled we can treat the motion of the molecule along the x, y and 
z directions independently. This justifies the use of the one dimensional model for the 
motion of the molecule. 

While for all cooling schemes in the optical domain the temperature of the 
electromagnetic environment can safely be set to zero, this assumption is in general 
not valid in the microwave regime. In a cryogenic environment with a temperature T ~ 
0.02— 4 K the thermal occupation number of the cavity mode, N = l/[exp(hu> c / k-BT) — l], 
is typically of the order of N ~ 0—10. In addition to the vacuum contribution in equation 
(J3J), a finite N results in cavity-assisted decay of rate 7(iV + 1) and a cavity-assisted 
excitation rate jN. As we discuss in more detail below, thermally activated processes 
are one of the limiting factors in achieving cold temperature. 

2.2. Cavity cooling of polar molecules in the Lamb-Dicke regime 

When trapped close to the surface of the chip, cooling forces on the molecule arise 
from variations of the microwave fields on the typical length scale d ~ 1 /im, which 
is only determined by the electrode geometry. For a harmonic trapping potential we 
compare this scale with the harmonic oscillator length xq = ^Jh/2mv and introduce 
the dimensionless parameters r\ = xoVQ/Q\ x= o and r) c = x V g / g\ x =o in analogy to the 
Lamb-Dicke parameter in optical cooling schemes; e.g. in the case of trapped ions, 
7] = 2nxo/X with A the laser wavelength [36j. As the linear extension of the molecule 
trap is approximately limited by d we find that cooling in this setup naturally occurs 
in the so-called Lamb-Dicke limit rj, r) c ~ Xo/d <C 1, where the spatial extension of the 
molecule wavepacket is small compared to variations of the microwave fields. While 
typically r\ and t] c will be equal or at least of the same order, in the following discussion 
we treat them as two independent parameters. 



h{A c - A)c*c, 




(4) 
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In the Lamb-Dicke limit we can expand the Hamiltonian (j3J) to first order in x, i.e. 
H$ ~ H + Hi. Here H is the Hamiltonian describing the free evolution of the internal 
and motional degrees of freedom of the molecule. The coupling between the internal 
and motional degrees of freedom is of the form 

H x = xF, F — + HVg(a + c + , (5) 

such that (F) is the mean dipole force exerted on the molecule. The small parameters 
r] and r] c allow us to treat Hi as a weak perturbation for the decoupled dynamics of 
the system. For |ifi| <C Ti, with T\ = min{7, k] the characteristic relaxation rate 
of the cavity-molecule system, we can adiabatically eliminate the internal degrees of 
freedom (see section [3]) and derive an effective equation of motion for the reduced 
external density operator fi(t). The resulting equation for the mean occupation number 
(n)(t) = Tr{aV0)} (see also [2H1 ESI S3) is of the form 

(n) = -W \{n) - (n) ] . (6) 

For a positive damping rate W > 0, equation ([6]) describes relaxation of the mean 
occupation number to a final steady state value (n)o- The relevant quantities 
characterizing the cooling dynamics, W and (n) , are given by 

W = A_ - A+, (n) = A + /(A_ - A+) , (7) 

where the cooling and heating rates A_ = S(u) and A + = S(—u) are determined by the 
spectrum of the force operator, 

S(u) = Re dr (F(r)F(0)) e + ^ . (8) 

Equations (ED-dHJ) provide the general relation between the cooling dynamics of the 
external degrees of freedom and the fluctuation spectrum of the internal operator F. 
Efficient cooling schemes correspond to the case where both the difference S(v) — S(— v) 
is maximized (large cooling rate) and the ratio S(—v)/S(v) is minimized (low final 
temperatures). Throughout this paper we make extensive use of this relation to discuss 
different cooling schemes in terms of the properties of S(uj). 



2. 3. Cavity cooling of polar molecules in the bad cavity regime 

While the results for cooling in the Lamb-Dicke limit summarized in equations (jH])-(jBD 
are valid for arbitrary cavity-molecule interactions, we concentrate in the remainder of 
this paper on the bad cavity limit (BCL), g k, which we expect to be the most 
relevant regime in experiments. For additional discussions on cavity cooling of atoms 
in the good cavity limit the reader is referred to [29J. The BCL allows us to eliminate 
the cavity degrees of freedom and to study and interpret cavity-cooling in terms of the 
resulting effective dynamics for the molecule. However, based on the underlying coherent 
interaction and the finite bandwidth of the cavity mode, the effective dynamics of the 
molecule differs significantly from the usual case of a two level atom coupled to the 
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Figure 2. Example for the force spectrum S(u>) for zero temperature (N = 0) and 
for a) ft <C |A| and b) Q ~ |A|. In both figures the contribution of Sq(uj) and S g (u>) 
are indicated by the dashed and dotted lines, respectively. 

free radiation field. The adiabatic elimination of the cavity mode, which is discussed in 
more detail in section [3], results in a decomposition of the general form of S(u) given in 
equation (jSJ) into three contributions, 

S{u) = S a (u)+S g {u) + S I {u). (9) 

The labels on the first two terms refer to their origin; the spectrum Sq is associated with 
the force from the drive field gradient Vf2 and exhibits sharp resonances of width ~ 7, 
whereas S g is associated with the force due to the cavity-TLS coupling gradient, Vg, and 
exhibits a broader resonance of width 2k. As explained in more detail in the following 
discussion, we can with each of these forces associate a different cooling mechanism 
in our system. The third contribution to the spectrum, Si, describes destructive and 
constructive interference between the different cooling processes. For illustration, the 
spectrum S(u) and its decomposition according to equation (jUJ) are shown in figure [2] 
for the case of zero temperature and a) for a weakly (Q -C |A|) and b) for a strongly 
(f2 ~ |A|) driven molecule. 

The results presented in this overview section apply to a weakly driven molecule, 
and generalizations to arbitrary Q are derived in section HI 

2.3.1. Cavity- assisted sideband/ Doppler cooling. We first consider the limit r] c — ► 
where the cavity field exerts no force on the molecule but still provides an effective 
decay channel for the excited state |e). In this limit the force spectrum reduces to 
the first term in equation (jHJ), S(u) = Sq(u). For zero temperature and fi C 7 the 
spectrum exhibits a single Lorentzian peak of width 7 centered at a frequency u = —A 
(see figure [3]). We identify Sn(u) with the excitation spectrum of a decaying TLS and 
recover a situation which is familiar from laser cooling of trapped ions [36l [37] . In the 
so-called resolved sideband limit, i.e. when the trap frequency is large compared to the 
effective decay rate 7, we can choose the detuning A = — v such that S{y) 3> S(—u). 
Under those conditions sideband cooling leads to final occupation numbers close to the 
quantum ground state, (n)o = (7/4z/) 2 . 
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Figure 3. Cavity-assisted sideband cooling, a) Spectrum Sq(u>) in the weak driving 
regime ft <C |A| for N = (solid line) and N = 0.5 (dashed line). The parameters 
chosen for this plot are (all in units of v): £1 = 0.1, A = — 1, A c = 0, K = 5, g = 0.5. b) 
Resonant processes contributing in the sideband resolved regime. Dashed lines indicate 
processes at finite environment temperature. 



For non-zero temperature the cavity acts as a finite temperature reservoir for 
the molecule, which has a corresponding equilibrium excited state population p° ee = 
N/(2N + 1). In the spectrum another resonance appears at oj — —v with a height 
~ Pg e . As seen from figure [3] b) this resonance corresponds to transitions from the 
excited state to the ground state which for A = — v preferably increase the motional 
quantum number and lead to heating. In equilibrium heating and cooling rates balance 
each other, meaning that the steady state occupation number (n)o in the trap is equal 
to the thermal cavity occupation number N. This equilibrium condition translates into 
a final temperature Tf = (u/uj c )T <C T. In section H] we show that for a weak driving 
field Q <^ 7jv = (2N + 1)7 and for 7^ <^ v the correct interpolation between the 
zero temperature limit (n) = (7/41/) 2 and the finite temperature result (n) ~ of 
sideband cooling is given by 

(n) = iV + (2iV+l)(^) 2 . (10) 

While in the sideband resolved regime cooling is most efficient we may instead consider 
the so-called Doppler limit 7^ 3> v. In this regime cooling rates are optimized by 
choosing the detuning A = — 7^/2 and the minimum final occupation number is, 

<n)„ = (2JV+l)g. (11) 

In contrast to the sideband resolved regime, we find that when the linewidth 77V exceeds 
the trapping frequency v the scaling of the final occupation number with temperature 
is less favorable (~ A^ 2 ) and cooling becomes inefficient for large N. 

Equations fllOp and (ITTj) have been derived for a weakly driven molecule and 
generalize the usual limits for sideband and Doppler cooling of trapped ions to the 
case of a finite temperature reservoir. In section H] we study the limit of a strongly 
driven molecule (see figure [2] b)), where we show that by increasing Q we benefit not 
only from a higher damping rate W, but cavity assisted sideband cooling also gets more 
robust with respect to imperfections in the trapping potential. 



Theory of cavity-assisted microwave cooling of polar molecules. 



10 




Figure 4. Spectrum S g (co) in arbitrary units for the set of parameters (all in units of 
v): g = 0.01, k = 0.3, fl = 2, A = 6, A c = 7, N = 0.5 b) processes which contribute 
to the individual resonances in the sideband resolved regime. 



2.3.2. 'Vg-cooling. We now consider the opposite limit r/ — ► where cooling forces 
arise from gradients of the cavity-molecule coupling g(x) only. In addition, to avoid 
interference effects discussed below, for the moment we assume g — > such that the 
spectrum offeree fluctuations reduces to the second term in equation Q, S(uj) = S 9 (uj). 

For simplicity we first look at the case of zero temperature N = 0. As shown in 
figure El S g (uj) has a single resonance at uj = A c — A with a width equal to the photon 
loss rate 2k. The deviation from the naively expected value uj = A c can be understood 
from the level diagram shown in figure H]b). With the molecule initially in the ground 
state, the energy hu extracted from the trap has to match the energy to fully excite 
the TLS (—HA) and the additional energy to excite the cavity (hA c ). In analogy to the 
discussion given above we distinguish two regimes of Vg-cooling, namely the resolved 
sideband regime k < (SB) and the Doppler regime k ^> v (D) with the corresponding 
limits for the steady state occupation number, 

(SB): {n)*=(£)\ (D): <n)„ = (12) 

Therefore, at zero temperature the final occupation number produced by Vg-cooling 
may be discussed in analogy to the cavity-assisted cooling, with the decay rate 7 replaced 
by the photon decay rate 2k. The sideband resolved regime, which has been discussed 
in the optical regime e.g. in reference [H], seems to be accessible in the present setup. 

For finite temperature, Vg-cooling becomes more involved as can be seen from 
the spectrum S g (u) plotted in figure HI In the resolved sideband regime, S g (u)) now 
consists of four distinct resonances at u = ±(A C — A) and u = ±A C . The resonance 
at w = A c — A corresponds to the process discussed above for zero temperature, 
and the resonance at uj = — (A c — A) corresponds to the analogue process with the 
molecule initially in the excited state. The two resonances at the cavity detuning 
uj = ±A C describe thermally activated processes which do not involve the external 
driving field Q. We interpret the corresponding heating and cooling rates as a cavity 
mediated thermalization process of the molecular motion, which otherwise would be 
highly decoupled from the electromagnetic environment. In the weak driving limit we 
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find that the resonances at u = ±A C are of equal height and therefore result in a purely 
diffusive (W — 0) dynamics. In other words, due to the large frequency mismatch 
oj c ^> v implicitly assumed in the derivation of our model, any finite thermal occupation 
of the cavity mode translates into an effective infinite temperature reservoir for the 
molecular motion. 

In the resolved sideband limit k C y we minimize re-thermalization rates by 
choosing either |A C | ^> |A C — A| or |A C | = 0. In section H] we show that Vg-cooling at 
finite temperatures results in a steady state occupation number, 



where (ra)th ~ 8(iV + 1) p° ee {i% / Vt) 2 for both choices on |A C |. Thus in contrast to cavity- 
assisted sideband cooling, we find that for Vg-cooling at finite temperature the final 
occupation number depends crucially on the driving field Q. A similar result is obtained 
for the Doppler regime, 



with /9 th ~ 16 (N + l)pg e (fi;/f2) 2 for the same choices on |A C |. From equations (fT3j) and 
(fT4l) we conclude that to employ the mechanism of Vg-cooling efficiently even at finite 
temperature we must go beyond the weak driving limit Q — > and consider either 
~ k or Q ~ |A|. In the regime fi~K<C|A|,|A c |, the cavity is far detuned from the 
rotational transition frequency and cavity mediated thermalization is suppressed, but it 
requires enough power to scatter photons from the driving field off-resonantly into the 
cavity. Alternatively we may consider the regime Q ~ |A| such that the driven excited 
state population p^ e = fi 2 /[(4A 2 + 2Vt 2 + 7 2 )] becomes comparable to p° ee = N/(2N+1) 
and stimulated cooling transitions eventually dominate over thermally activated heating 
processes. However, in this regime equation (1131) and ( 1T4"1) are no longer valid and we 
postpone the discussion of Vg cooling in the strongly driven regime fl |A| to section 
H] where we find a significant improvement also for a low driving power Q > 7. 

2.3.3. Cooling by quantum interference. As a consequence of interference effects 
between forces exerted on the molecule by the cavity and the classical microwave field, 
the spectrum S(u) plotted in figure [2] clearly differs from the naively expected sum of 
Sq(u) and S g (u>). In equation ([9]) we have introduced the term Si(uj) to account for 
this difference. 

In figure [5] we have sketched the three different paths which contribute to the 
interference effects. In a) the molecule is excited on the red sideband followed by an 
effective decay via the cavity on the carrier transition. This process is described by 
Sq(uj). In b) the molecule is excited on the carrier transition and changes its motional 
state in the decay process, described by S g (u). The third path sketched in figure [5] 
c) might not be obvious in the first place and has the following interpretation. The 
position dependent coupling of the TLS to the cavity mode, g(x), in general leads to a 
position dependent decay constant j(x) as well as a position dependent Stark shift S(x). 




(13) 




[(2A+l)+/3 t h] 



(14) 
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Figure 5. a)-c) Different paths which interfere in the cooling process. To emphasize 
the underlying coherent interaction we have decomposed the effective decay process 
in a coherent energy transfer (g) and a successive emission of a cavity photon, d) 
Contribution to the cooling rate which at finite temperature adds incoherently to the 
other paths. 



In a semiclassical picture, the oscillating motion of the molecule in the trap therefore 
translates into modulations of 7 and 5 at the frequency v and results in additional 
sidebands in the excitation process. Note that interference between path b) and path 
c) also occurs for 77 — > and in general we cannot write S(u) = S 9 (lj) even in this 
limit. However, in section H] we will argue that in the parameter regimes relevant for 
Vg cooling those interference effects only slightly modify the cooling process. For finite 
temperature a fourth process, the decay of a thermally excited molecule as sketched in 
figure [5] d), adds incoherently to the processes a)-c). Again, this process is proportional 
to the thermal equilibrium excited state population p° ee . 

Cooling by quantum interference was discussed and demonstrated in the context 
of electromagnetically induced transparency [IS, H7] and has been proposed for atoms 
in optical cavities [261 12H]- In all those schemes cooling relies on a reduction or ideally 
a cancelation of blue sideband transitions (S(—u) — > 0) by interference, while S(u) 
is still finite, which in the ideal case leads to ground state cooling, (n) = 0. From 
figure [2] we see that in the present system such an exact cancelation of the heating 
rate appears in the weak driving regime, but is no longer present in the strongly driven 
system. At finite temperature thermally activated processes reduce the contrast of the 
interference features, and from the discussion of the individual paths given above we 
conclude that cooling by quantum interference is only efficient as long as driven cooling 
and heating processes ~ p^ e dominate over thermally activated processes ~ p° ee . Due 
to the sensitivity to temperature as well as other system parameters, we do not discuss 
cooling by quantum interference in much detail in the present work, but in section H] 
we point out several consequences of interference effects for the combination of CASC 
and Vg-cooling. In references [26| [28l [29] more details on interference effects in cavity- 
assisted cooling of atoms can be found for the weak excitation regime (N — 0, Q — > 0). 
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3. Theory of Cavity Cooling in the bad cavity limit 

In the previous overview section we have presented for the case of a weakly driven 
molecule a physical picture of the individual cooling mechanisms which contribute to 
the overall cavity-assisted cooling process. For the weak driving regime and for zero 
temperature, (f2 — »• 0, N = 0) analytic expressions for S(u) and the resulting cooling 
rates have been derived in previous work [261 I2H] on a related system of trapped atoms 
inside optical cavities. However, for non-zero temperature we have identified thermal 
cavity-mediated heating processes which degrade the cooling efficiency, and thus in the 
present setup there is a clear necessity to extend the discussion of cavity cooling schemes 
to the regime of a strongly driven molecule. In this section we present a derivation of 
heating and cooling rates which is valid in the bad cavity limit (j<k and allows us to 
study cavity cooling for a wide range of system parameter. 

The starting point for the derivation is the combined cavity-molecule system 
discussed in section [2} for which the evolution of the total density operator p(t) is 
determined by the following master equation, 

p = C(p) = -i[H s ,p] + C K (p). (15) 

The system Hamiltonian H$ is defined in equations ([I]) and (j3J), and the coupling of the 
cavity field to a thermal reservoir is taken into account by the term 

C K (p) = k{N + 1) V[c}(p) + kN Z>[ c t](p) , (16) 

with V[c](p) = 2cpc { ' — pc^ c— cr cp . Starting from the full master equation (115j) we proceed 
in two steps. First, we make use of the small parameter g/n to eliminate the dynamics 
of the cavity mode and derive an effective master equation for the internal (TLS) and 
external degrees of freedom of the molecule. In a second step, the weak coupling between 
the TLS and the external degrees of freedom, represented by the small parameters rj, rj c , 
also allow us to eliminate the TLS and derive an effective master equation describing 
only the motion of the molecule in the trap. 

3.1. Adiabatic elimination of the cavity mode 

For the elimination of the cavity degrees of freedom we decompose the total Liouville 
operator £ defined in equation (flET) into three contributions 

C = C c + C m + C g . (17) 
The first two terms describe the uncoupled dynamics of the cavity mode, 

C c (p) = -i[{A c - A)c t c, p] + C K (p), (18) 
and the molecular degrees of freedom, 

C m {p) = -i[H mi p], (19) 

whereas the interaction between the cavity mode and the molecule is described by the 
last term in equation (fT7|) . 

C g (j>) = -i\g(x)(<T+c+<T-J),p]. (20) 
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In the decoupled limit g — > the cavity field evolves independently of the molecule. 
On a timescale k" 1 the system relaxes into the state p(t) ~ p° c ® p m {t), with p° c the 
cavity equilibrium density operator defined by £ c (p°) = 0, and the molecular operator 
Pm.it) evolving under the action of £ m . For a finite coupling j < k deviations from 
the factorized form of pit) are small, but the coupling term C g modifies the dynamics 
of Pm(t)- To proceed we adopt a projection operator technique along the lines of 
reference jl5] and define the projector Vp = p® <8> Tr c {p} and its orthogonal complement 
Q — (1 — V). Inserting the decomposition p(t) = Vp(t) + Qp(t) into equation (fT5"j) we 
obtain the two coupled equations 

Vp(t) = VC m Vp{t) + V£ g Qp(t), (21) 
Qp(t) = Q(C C + £ m + £ g )Qp(t) + QC g Vp(t). (22) 

As the population in the subspace Qp is damped with a rate k which is fast compared 
to the coupling term QC g V ~ g, we formally integrate equation (1221) . insert the result 
into equation ( T2TT) and expand the final expression up to second order in g. For times 
t ^> ft -1 we end up with an effective master equation for the molecule density operator 
Pm(t) = Tr c {V p(t)} which is given by 

POO 

p m (t) = -i[H m ,p m (t)} + / drTi c {C g Qe^ +c ^ T QC g e- c - T ip m {t) <g> p° c )}. (23) 

Jo 

Strictly speaking the effective master equation (123]) is valid in the limit g^/N + 1/k 1, 
but we expect it to hold to a good approximation also for a slightly relaxed condition, 
9 ^ «■ 

The second term in equation ( |23|) describes the effect of the cavity on the molecular 
dynamics. We evaluate this term by inserting the definitions of C c , £ m and C g given in 
equations (fTBl) - (l2"01) and obtain an effective master equation on the form (see also [43] ) 

p m = -i[H m , Pm ] + [(N + 1) (Tp m Si - S^T Pm ) + N (T^ Pm S - ST^p m ) + H.c] . (24) 

Here we introduced the abbreviations S = g{x)a- and 

POO 

T= dre- !(Ac - A)T e' KT S(-r), S{t) ^ e iHmt Se~ iHmt . (25) 
Jo 

In summary we have derived an effective master equation (1241) for the molecule, where 
the first term describes the bare evolution, and the second term describes the influence 
of the cavity on the dynamics. We want to point out that the latter term in general has a 
complicated form, since the operator T (1231) depends on the bare molecule Hamiltonian 
H m . Only in the infinite bandwidth limit, where k is large compared to A, Q and u, will 
equation (1231) reduce to the simple form T = S/(k — iA c ). 

3.2. Effective molecule dynamics in the Lamb-Dicke regime 

To proceed with the derivation, we now consider the master equation ( 1241) in the Lamb- 
Dicke limit r), rj c 1 where internal and external degrees of freedom interact weakly. 



Theory of cavity-assisted microwave cooling of polar molecules. 



15 



We expand the position dependent Rabi frequencies Q(x) and g(x) up to first order in 
the Lamb-Dicke parameters, i.e., 

Q(x) ~ +r/(a + a f )] and S ~ ga-[l + rj c (a + a f )]. (26) 

Under the same approximation the operator T gets the form 

T — g [£_ + r/ c S_(z/)a + 77 c E_(-i/)a+] , (27) 

Sl(— i/) effective TLS operators defined in (1A.3I) . Up to second 



V 



with £_ 

order in the Lamb-Dicke parameters we obtain the master equation 

Pm (t) ~ (C + Ci+ C 2 ) p m {t) , 



(2* 



where the subscripts indicate the respective order in the Lamb-Dicke parameters rj, r] c . 
The zeroth order term is given by 

Co = & + £ E , (29) 

where £e(p) = — iv[cfia, p] describes the uncoupled external dynamics of the molecule, 
and C\ describes the effective dynamics of the uncoupled TLS, 



-i[Hi,Pv 



+ g 2 N 



-PmCr 



S_p m CT + 



cr + E_p E 



H.c. 



H.c. 



(30) 



Explicit expressions for C\ and £2 are listed in Appendix A Let us for the moment 
consider the effective dynamics of the TLS ( l30i) . with the corresponding Bloch equations 
for the set of operators {o~ x , a y , a z } of the form 



(31) 
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Here 7^ = (2A^ + 1)7 and the remaining parameters are defined in Appendix B In the 
discussion in section H] we will mainly focus on the limit Q k where the parameters 
S x ,y, Jx,y, &>x,y and T x<y are small and equation ( 131]) reduces to the normal Bloch equations 
for a two level system coupled to a thermal bath. Further the stark shift 5 and the 
effective decay rate 7 reduce to the familiar expressions 



8 -> 5 = (2N+ 1) 



g 2 \ 



K 2 



A? 



7^7 



2p 2 K 
■? + Ai 



(32) 



For the purpose of readability, we will use the simplification A + 5 — ► A in the following. 
Finally note that when the driving strength Q is comparable to the cavity linewidth 
k, unless |A|, |A C | ^> k there will be modifications of 5 and 7. In addition, the non- 
standard terms in the Bloch equations ( 131]) introduce new features in the TLS dynamics 
which modify the force spectrum S(uu) and in the end the cooling dynamics. A detailed 
discussion of the effective dynamics of a strongly driven TLS inside a cavity is given in 
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3.3. Adiabatic elimination of the internal degrees of freedom 



Starting from equation ( 128]) we further eliminate the internal degrees of freedom in a 
similar fashion as the above elimination of the cavity mode, following the procedure in 
[35] . First we notice that the zeroth order term £ in the expansion (1281) describes the 
uncoupled dynamics of the internal and external degrees of freedom, whereas the higher 
order terms C,x t i couple the external and internal dynamics and therefore introduce 
forces acting on the molecule motion. Provided the coupling is weak compared to the 
energy scales defining the TLS dynamics, r]Q,rj c g 2 /K, <C 7tv, the state of the TLS will 
deviate only slightly from the equilibrium state due to the interaction, and the dynamics 
of the TLS can be eliminated. We define a projection operator V^, 

^(Pm) = Pi ® ^ l 71 ^ 71 ! Tr l{H/°m|™)} , (33) 

n 

with pj the internal equilibrium state defined by £o(p?) = 0- Second order perturbation 
theory in the Lamb-Dicke parameters t], 7] c <C 1 results in a master equation for the 
external operator fi(t) = Tr I {'P At p m (t)}, which describes the evolution of the population 
of the individual trap levels. We obtain 

POO 

fi{t) = / dr^iVC^QCM®^))} + ^iOV^Pi ® M*))}- (34) 
Jo 



Inserting the definitions of £i ; 2 given in |Appendix A| we end up with a master equation 
on the form 

fi{t) = AJD[a] {pit)) + A + V[a)} (//(*)) , (35) 

and recover the cooling equation ([6]) given in section [2] with W = (A_ — A + ) and 
(n)o = A+/W. The heating and cooling coefficients A± which follow from an evaluation 
of equation (]34p can be divided into three contributions, 

A ± = S n {Tv) + S 9 {tv) + Sj(Tiy) , (36) 

where the physical motivation for the decomposition ( 1361) has been presented in section 
[2j The first contribution to the force spectrum, Sq(uj), follows from the first term of 
equation (1341) in the limit t] c — > and is defined as 

Sn(u)) = ■!— Re / drTiiKe^M)}^ . (37) 
* Jo 

The second term, S g (u), follows from the second term in equation (1341) and is given by 

S g (u) = 2rf c g 2 Re{(iV + I)(t7 + E_(w))i + N^JZ+lw))^, (38) 

with effective TLS operators Yi±(y) given by flA.5[) . Finally, remaining contributions 
from the first term in equation ([Ml) are found by making the Ansatz, 

®fi) = ® W + £-(Pi) ® «V + ^L(Pi) ® P« + £+{Pi) ® A*a f • (39) 

The result is summarized by Sj(a>) defined as 

/■oo 

Sj(±u) = 2Re / rfrTr I {/Cl(e £( ' r /C ± (p°))}e ± ^ - S a (±u), (40) 
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with the action of the superoperators K,± given by 

K>±{p) = -i^°~xP + r) c g 2 {N + l)[cr_p£ + + £_(±^)P°~+ - o- + S_(±^)p - er + £_p] 

+ r] c g 2 N[a + p^ + £ + (±z/)pa_ - a„S + (±z/)p - <r_£ + p] . (41) 

In conclusion, we have expressed heating and cooling rates A± in terms of two-point 
correlation functions and steady state expectation values of a TLS whose evolution 
is described in terms of the effective Bloch equations ( 13TT) . The results are valid for 
gy/N + 1 <C k and rjVt <C jn, r] c <C 1 but at this stage no further assumptions have 
been made. However, in the following we mainly concentrate on the case O < k or on 
the regime ~ k, |A c |, |A| ^> Q where the Bloch equations reduce to the standard form 



,see 



Appendix B[). Using the quantum regression theorem and the steady state solution 



of equation ( 13TT) we outline the derivation of analytic expressions for Sn(u) and S g (u) 
Appendix C and Appendix D Analytic expressions for Si(u) could in principle be 



111 



derived along the same lines but in general do not have simple enough form to provide 
further insight. Therefore we will rather use numerical calculations in combination with 
previously derived results for the weak excitation limit [29] as a basis for the discussion 
of the interference terms. 



4. Discussion 



In this section we discuss cavity-assisted cooling of polar molecules both in the weak 
and strong driving regime, and base the discussion on the result for S(u) derived in 
section G2 For a physical understanding we again find it useful to discuss cavity assisted 
sideband cooling and Vg cooling independently in a first step. In a second step we will 
then study the general case where both mechanisms as well as interference effects are 
taken into account. Analytical results for cooling rates and final occupation numbers 
for the different parameter regimes are summarized at the end of this paper in table [TJ 



4-1. Cavity assisted sideband cooling 



We start the discussion with the mechanism of cavity assisted sideband cooling (CASC) 
where the corresponding cooling rate Wq, = Su(u) — Sq(— v) and steady state occupation 
{ n )o,n — Sq(— u)/Wq depend on the spectrum Sn(u) defined in equation (j37]) . Following 
the calculations outlined in Appendix C we obtain for the spectrum in the weak driving 
limit, Q 7at, 

V 2 Q 2 



+ 



(42) 



L(o; + A)2 + 7 y4 (u;-A) 2 + 7 y4_ 

For a fixed ratio v/^n we optimize (n) ^ with respect to the detuning A, where 
obviously A must be negative to promote cooling. For the two limits of interest, the 
resolved sideband limit (SB) v ^> 77V and the Doppler limit (D) v <ti 77V, we find A = —v 
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and A = —-f^/2 respectively. The final occupation number is given by the respective 
expressions ffTUl) and ffTTj) in section HJ and the corresponding cooling rates are given by, 

(SB) - Wq -(2N + 1) 1n > (D) - Wa ~ (2N+1) 1N - (43) 
Although in the weak driving limit the final occupation number is independent of Q, 
the cooling rate scales as Wq ~ fi 2 and under the validity of equation (H31) it is bounded 
be Wq < t] 2 '-/. In competition with the cavity-mediated thermalization process of the 
molecule, discussed in section [2} it is therefore necessary to study sideband cooling 
beyond the weak driving regime. 

Let us now consider an arbitrary Q, but focus on CASC in the sideband resolved 



limit 7^ < 1. As shown in Appendix C, the spectrum now exhibits a three-peak 
structure and reads, 

v 2 n 2 



a+7 "o7o «-7 



(44) 



^-A) 2 + 7 2 /4 w 2 + 7^/4 (cj + A) 2 + 7 2 /4_ 

with A = \/Vl 2 + A 2 . For sideband cooling we are only interested in the peaks at nonzero 
frequency. The relative heights of the resonances, a±, and the widths 7,70 which are 



proportional to 7^ are defined in Appendix C and are functions of f2, A and N. For Wq, 



to be positive we require a red detuned microwave field A < and cooling is optimized 
for v = A. Under those conditions we obtain 

^.s^^j.^y^, (45) 

with g( if) = 4 sin 2 tp | cos 3 cp\ /(4 — sin 4 ip) , and 

(n>o,n = ^ + 0(1n/" 2 ) * N + (2N + l) gl (^, (46) 

withgi(o?) = (1 — I cos<^|) 2 /(4| coso?!). The functions g(<p) andgi(y9) are plotted in figure 
E]and depend on the angle ip defined as sin(o9) = Q/A. We find that the cooling rate 
increases with Q up to Q ~ |A| but then goes to zero for Q 3> |A|. It has a maximum 
value of g((po) — 0.2 at smip = Qq/A ~ 0.65, and we want to point out that this 
holds independently of the temperature N. The final steady state occupation number 
increases monotonically with increasing Q; gi(v?o) — 0.02 and gi — > 00 for Q ^> |A|. 

We conclude that driving moderately, Q ~ |A| ~ u/2, increases the cooling 
rate tremendously while only marginally increasing the final occupation number. We 
point out that the validity of the adiabatic elimination constraints the parameters, 
7]v ~ r]Q < 7tv, and thus the expected upper limit W < 7 still holds. Nevertheless, this 
bound is by a factor I/77 2 higher than in the weak driving case. 

4-2. Vg cooling 

Let us now move on to discuss the mechanism of Vg cooling, where the cooling rate 
W g = Sg^v) — S g (—v) and t\iQ final occupation number {n)o, g = S g {—u)/W g are given 



in terms of the spectrum S g {u) defined in equation (|38|) . In Appendix D we show that 
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Figure 6. Dimensionless cooling rate g(<p) and final occupation number gi(<p) for 
CASC in the resolved sideband limit, as functions of sintp = VI/ v. The vertical line 
marks the optimal drive strength VI = 0.65^ which maximizes the cooling rate. 



for finite temperature N ^ and arbitrary f2 k or Q ~ |A|,|A c | the structure of 
S g (u) flD.41) is quite involved and generally exhibits 6 resonances centered at frequencies 
oj = ±(A C — A) and uj = ±(A C — A ± A). Note that we assume 7 at <C A, which holds 
naturally when we later combine Vg cooling with CASC. As the height of two of the 
resonances, at uj = ±(A C — A — sA) where s = sgn(A), scales as ~ (f2/A) 4 we can for 
weak driving Q <C |A| neglect those contributions and obtain the spectrum 

c L ,-^Vr ^ + i)(W kN(q/a) 2 v MN + l)p° ee 1 , , 

g[ ' ~ 2 I k 2 + (u- A g ) 2 k 2 + (uj + A g ) 2 ^k 2 + (uj + sA c ) 2 \ " 1 ' 



We note that only the driven terms in equation fl47l) contribute to the cooling rate W g , 
which for optimized detuning A g = A c — A = n in the Doppler limit and A c — A = v 
in the sideband limit, respectively, are given by 

(SB): g)* (D): (48, 

The last term in (1471) is symmetric in frequency and thus corresponds to a purely 
diffusive process which is present only when the thermal excited state population p° ee is 
nonzero. For finite temperature there is thus a competition between the driven cooling 
processes proportional to p ee ~ Q 2 /A 2 and thermal diffusive processes. In the resulting 
expressions for the final occupation number (n) 0t9 , which are given in equations (jl~3r) 
and (jl~4]) . corrections due to the additional diffusion process appear in the quantities Ah 
and (n)th which in general are given by 

(n) - 4^+lV ( {A/Q)2 + (A/ ^H 3 - 1 ^ 2 ( Ar + 1 )^(A/^) 2 , . 
(n) th - 4k (N+l) Pee ^ 2 + (Ac + A)2 + ), Ah - -TT^ A 49 ) 

As already mentioned in section El to reduce Ah and (n)th there are two different 
strategies. We can either stay in the far detuned regime | A c | ~ | A| ^> k, Q and increase 
the driving power Q ~ k, or alternatively, use a resonant driving field |A| — > 0. In the 
first case we are still in the far detuned regime, Q -C |A|, and equation ( 14TI) is still valid, 
even for n ~ Q. For a resonantly driven molecule we however need to reevaluate S 9 (uj), 
as we will do in the following. 
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We now consider the strong driving regime Q ^> |A|,7jv where the rotational 
transition of the molecule is saturated; p ee — ► 1/2. In this limit the general spectrum 
(1D.4I) is of the form, 

2k(N + 1) 



Sg(uj) 



Vc9 2 
4 



+ 



Vc9 2 
4 



«2 + [u> - (A c - A)] 2 
2kN 

k 2 +[u + (A c - A)] 2 



E 



=±i 



w(jV+lj 
[w - (A c H 

kN 



s=±l 



K 



(A c + sQ)]< 



(50) 



We find that for a saturated molecule at moderate driving strength |A|,7jv < < 
where the second condition is already implicitly assumed in the derivation of equation 
( l5T)j) . the three peaks overlap almost completely and thus to a very good approximation 
we can consider all resonances to be centered around A c . We then obtain the two peak 
structure 

k(N+1) kN 



2 2 
Vc 9 



+ 



(51) 



K 2 + (iO-A c ) 2 ' k 2 + (uj + A c 

in analogy to the spectrum Sq(u) (142 jl in the weakly driven regime. In the fully saturated 
regime the final temperatures achievable by Vg cooling are given by 



(SB) : (n) 



0,9 



N 



( 2N + 1 ^ 



(D) : (n) , g 



(2iV+l)f , 
2v 



(52) 



and exhibit a similar dependence as in the case of CASC but with 7 replaced by 2k. 
However, for a finite environment temperature (n) 0tg grows in both regimes only linearly 
with N. The cooling rates are independent of N and approach the optimal values of 

V 2 c 9 2 (v/k) 



(SB) : W g 



Vc9 2 



(D) : W g 



(53) 



K K 

In conclusion we find that the apparent inefficiency of Vg-cooling at finite 
temperatures can be overcome either in the regime where the cavity is far detuned 
from the rotational transition and the driving strength Q is comparable to the cavity 
linewidth k, or in the fully saturated regime A = 0, Q ^> 7. In the later case the cooling 
rates are optimized and the resulting final temperatures grow only linearly with N. 



4-3. Combining cavity-assisted microwave cooling and Wg cooling. 

Until now we have only considered the idealized limits 77 — ► and t] c — >• where either 
the mechanism of cavity- assisted sideband/Doppler cooling or Vg-cooling dominate the 
cooling process. In this section we focus on the more realistic case where both Lamb- 
Dicke parameters are of the same order and discuss the combination of the cooling 
mechanisms, being aware that in general W 7^ Wq + W g due to interference effects. In 
view of the tight confinement required to trap the molecule close to the surface of the 
chip, the following discussion is restricted to the resolved sideband regime 7^ < v. 

From our analysis so far, we conclude that under ideal conditions cavity cooling of 
polar molecules is dominated by the mechanism of CASC, which provides the highest 
cooling rates and results in the lowest final temperature both in the weak and strong 
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driving regime. However, although CASC is the more effective of the two mechanisms, 
we expect it also to be more sensitive to imperfections in the system, which in a chip 
based setup might be harder to avoid than in e.g. a large ion trap. Therefore, in 
the following we are particularly interested in the robustness of cavity cooling under 
non-ideal trap conditions. In order to discuss different types of imperfections like 
anharmonicity or a state dependent trapping potential within a general and simple 
model, we proceed as follows. We assume that the system parameters are optimized 
with respect to the expected trapping frequency v which differs from the real trapping 
frequency v r = v + 8v by a shift 5u, and then we study the dependence of the 
cooling parameters W(5v) and For an anharmonic trap we can then interpret 

W n = W(5v n ) as a level dependent cooling rate, assuming a frequency spacing of the trap 
levels u n = v + 8v n with 5uq = 0. Small deviations from a state-independent trapping 
potential (e.g. assuming v is the trap level spacing as seen by state \g) and v + 5u 
as seen by state |e)) translate into a level dependent detuning A n — > A — nSu. Since 
for sideband cooling, which is most sensitive to imperfections, a shift of the transition 
frequency and a shift of the trapping frequency affect cooling the same way we can 
interpret W n = W(nSu) as a level-dependent cooling rate as well. 

4-3.1. Weak saturation regime. We first consider the weak saturation regime Q <C |A| 
and assume the relation 7 < v < k for the damping parameters. We choose A = — v to 
maximize excitations of the molecule on the red sideband and choose A c = k—v with the 
intention to optimize Vg-cooling in the Doppler limit. In that case the spectrum S(ui) 
is similar to the one plotted in figure [2] a) in section [2j As an example we plot in figure 
[7] the cooling rate W(5u) and the final occupation number {n)o(5v) for a specific set of 
parameters. We find that in the weak saturation regime the cooling rate is dominated 
by the mechanism of CASC and under non-ideal conditions 5u ^ it scales as 



The correction due to the gradients of g(x) lead to an increase of W(5u) for 8v < 0, 
and a decrease for bv > 0. The asymmetry shows that corrections arise from 
interference effects in the excitation of the molecule. As expected from the small 
value of ~ Q 2 /A 2 <C 1, the pure mechanism of Vg-cooling plays essentially no 
role for the cooling rate of a weakly driven molecule. However, as discussed above, the 
mechanical coupling to the cavity introduces an additional source of diffusion which 
becomes important when the thermally activated excited state population p° ee exceeds 
Pg e - I n contrast to pure CASC, where {n)o(Su) is highly insensitive to 5u (since the 
heating rate scales like the cooling rate, A+(8v) ~ 1/[1 + {25u/^n) 2 ]), the additional 
diffusion for t] c > scales differently with 5v and will affect the final occupation number. 
In the limit k ^> v we obtain 




(54) 




(55) 
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Figure 7. Cavity cooling in the weak saturation regime for the parameters (all in 
units of v) Q, = 0.1, A = —1, k = 4, A c = 3, g — 0.2 and r\ = rj c . a) Cooling 
rate W(5v) in units of Wo = 7/ 2 f2 2 /7. The dashed (dashed-dotted) line indicates the 
corresponding result for r/ c = (rj = 0). b) Steady state occupation number (n)o(Sv) 
for N=0.1 (blue), N=0.5 (red) and N=l (green). Dashed lines for i] c = 0. 



with 



0,9 



2 2 

Jn_ 
4z/ 2 



(2N + 1) + 



Q 2 



(56) 



In summary we find that for a weakly driven molecule the cavity-assisted microwave 
cooling is dominated by CASC and under non-ideal conditions \8v\ > 7at the total 
cooling rate decreases as W ~ Wq ~ {^n/8v) 2 . For finite temperature the additional 
diffusion caused by thermal Vg processes affects the final occupation number and 
(n)o ~ N is reached only for 5i/ < 7 and (7^ /fi) 2 p° ee < 1. 



4-3.2. Strong driving regime. We now consider the case of a strongly driven molecule 
Q ~ |A|, where according to our previous discussion the cooling rate for CASC, Wq, is 
(close to) optimized. In figure [8] a) we plot the cooling rate W(5v) with a fixed A = v 
(A < 0) but different values of Q. Apart from the improvement compared to weaker 
drive under ideal conditions, we note that W is also less sensitive to deviations 8v for 
Q ~ |A|. We attribute this feature partially to power broadening of the resonances, 
and also to a stronger contribution from interference effects. The asymmetry of the 
latter improves the robustness of cavity cooling for b~v < but has a negative effect 
for 8v > 0. In figure [H] b) we also plot the final steady state occupation number for 
N = 0.5 and different values of Q. We see that by increasing Q up to f2 the improved 
robustness of the cooling rate translates into a corresponding robustness of the final 
temperature with respect to deviations from ideal harmonic trapping conditions. In 
contrast to Wn, the thermal diffusion rates do not increase with Q, and thus their 
contribution to the final occupation number is small; the limit (n)o = N can be reached 
for ('J%/v 2 )(N+ l)p°e < 1- When we increase Q beyond the optimal value of flo ~ 0.65u 
both the efficiency and the robustness of cavity cooling decrease again. 
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Figure 8. Cavity cooling in the strong driving regime for the parameters (all in units 
of v) A = 1 (A < 0), k = 4, A c = 3, g = 0.2 and rj = tj c . a) Cooling rate W{8v) in 
units of Wo = V 2 1 f° r N = and different values of SI = O.i (solid), 0.65 (dashed) , 
0.95 (dotted), b) Steady state occupation number (n)o(Sv) for N = 0.5. 




Figure 9. Spectrum S(u>) for Vg cooling in the saturated regime. The three different 
lines show S(u>) in units of So = i] 2 g 2 j 'k for A = (solid), A = f2 (dashed) and 
A = — O (dotted). The remaining parameters for this plot have been chosen as (all 
in units of v): g = 0.3, k = 3, A c = 3, £1 = 0.2 and r\ = r/ c . Multiple resonance in 
the frequency range |w| < \/ A 2 + il 2 are due to excitations of dressed states of the 
strongly driven TLS. 



4-3.3. Fully saturated regime. We finally consider the fully saturated regime f2 3> 
7at, |A| where the effect of CASC vanishes while Vg cooling is most effective. When 
we assume in addition that A ~ Q u, excitations on the motional sidebands do not 
play a relevant role any longer and we recover the case of pure Vg cooling as discussed 
above in section 14.21 According to that discussion, we conclude that cavity cooling in 
this limit has the same efficiency and robustness as CASC in the weak driving limit, 
but with the relevant scale 7at replaced by 2k. As 7a? <C k we trade cooling efficiency 
for robustness by switching between CASC and Vg cooling. 

The arguments presented in the previous paragraph are valid under the assumption 
v 3> Q ^> 7tv, |A|, but these conditions might not always be satisfied in a real 
experiment. Especially, if the trapping potential is slightly state dependent, the detuning 
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Table 1. Limits for the final occupation number (n)o and the cooling rate W for the 
different parameter regimes studied in this paper. For \7g cooling in the weak driving 
regime, optimized values for (n) t h,Pth ~ N(n/Q) 2 are given in section[2l 



"CASC" 


sideband limit (y 3> Jn) 


Doppler limit {y <C jn) 


weak driving 
regime 


(n) = N+(2N + l)( lN /Av) 2 
W = V 2 n 2 /((2N+l) lN ) 


(n) = (2N + l)~/ N /4v 

W = 2? 1 2 n 2 v/((2N + l) 1 2 N ) 


strong driving 
regime 


{n^^N + O^l/u 2 ) 
W ~ 0.2(?7iV7jv) 2 7 




"Vg cooling" 


sideband limit (y 3> n) 


Doppler limit {y <C k) 


weak driving 
regime 


(n) = N+{N+l)( K /2 V ) 2 +{n) th 
W = V 2 g 2 /(2n)(n/A) 2 


(n) = [(2N+l)+pt h ](K/2v) 
W = r)*g 2 v/{2K 2 )(Sl/A) 2 


saturated 
regime 


{n) = N+(2N + l){K/2v) 2 
W = rj 2 g 2 /k 


(n) = (2JV+l)(/c/2i/) 
W = r] 2 c g 2 vlK 2 



A depends, in a semiclassical picture, on the position of the molecule in the trap. In 
figure[9]we plot the total force spectrum S(lj) for the case A = and for A = ±Q. Under 
resonance conditions (note that our definition of A includes the cavity induced Stark 
shift 5) the spectrum is to a good approximation given by S g (u), assuming 7^ v. For 
A comparable to fl we find that S(±v) is reduced since the excited state population 
drops below its maximal value of 1/2. In addition, for red detuning, A < 0, destructive 
interference for uj > A and constructive interference for uj < A degrade Vg cooling, 
while a blue detuned driving field, A > 0, has just the opposite effect. We conclude 
that cavity cooling in the saturated regime can be efficiently applied for a range of 
parameters fulfilling the condition v > Q > 7^, A > 0. Outside this parameter regime, 
excitations of the molecule on the red and blue motional sidebands as well as destructive 
interference effects play an important role and Vg cooling no longer provides a valid 
description of the full cooling (or heating) dynamics. 

5. Summary and conclusion 

In this paper we have studied cavity-assisted cooling methods for a single polar 
molecule, which are based on a strong coupling between two rotational states of the 
molecule and an on-chip microwave resonator. We have shown that in the Lamb-Dicke 
limit the cooling dynamics in this system can be understood in terms of two main 
cooling mechanisms, namely cavity-assisted sideband cooling (CASC) and Vg cooling, 
interference effects between those two mechanisms and a cavity mediated thermalization 
channel introducing additional heating at finite temperature. Our main conclusions are 
as follows. 

Cavity-cooling is most efficient by employing CASC, meaning that the molecule is 
resonantly excited on the red motional sideband followed by cavity-enhanced decay of 
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rate 7. In the sideband resolved regime, 77V = 7(2iV + 1) <C v, the lower bound for the 
final occupation number (n) is then given by N, the thermal equilibrium occupation 
number of the cavity mode, and translates into a final temperature of Tf = (u/uj c )T <C T. 
The cooling rate, which scales as ~ Q 2 for low excitation power is optimized for a value 
Q ~ 0.65^ where cooling rates close to the maximum value of W ~ 7 can be achieved. 
Including imperfections like the anharmonicity of the trapping potential characterized 
in our model by a level dependent shift of the trapping frequency Su n , the cooling 
rate decreases as W n ~ (7^ / '5u n ) 2 and eventually re-thermalization processes become 
dominant. Our analysis shows that deviations from the ideal case (n) — N are expected 
when either the frequency shift exceeds the excitation linewidth \b~v n \ > 7^ or when the 
sideband cooling rate drops below the re-thermalization rate for an insufficient driving 
strength Q. Together with the cooling rate also the robustness of CASC is optimized 
for Q ro O . 

In the far off- resonant regime | A|, | A c | ^> is, Q as well as under resonance conditions 
A pa cavity cooling is dominated by the mechanism of Vg-cooling. In the resonant 
case, Vg cooling is optimized in the fully saturated regime 3> 7jv, |A|, where the limit 
(n)o — N can be reached when the cavity linewidth 2n is smaller than the trapping 
frequency, while the result for Doppler cooling (n) = (2N + l)n/(2u) applies in the 
opposite limit. The maximal cooling rate under resonance conditions is W = r) 2 g 2 / k 
while in the off-resonant case the molecule is only virtually excited and the cooling rate 
is reduced by a factor 2 /A 2 . In this case, to overcome thermally activated heating 
processes the additional condition Q > k is required to reach a final occupation number 
close to N. While in general it might be harder to reach the sideband resolved limit 
k < ^, Vg cooling is by far less sensitive to imperfections as the condition 5u n < 7^ is 
replaced by b~v n < 2k. 

In this work we have analyzed cooling rates and final occupation numbers for a 
fixed set of parameters to obtain a physical understanding of the characteristics and 
limitations of cavity cooling in different parameter regimes. From this analysis more 
efficient cooling strategies can be derived by considering the possibilities to switch 
between different parameter regimes during the cooling process. For example, at the 
initial stage where the anharmonicity of the trapping potential may be most severe, Vg 
cooling in the Doppler regime offers a very robust cooling method. In a second stage, 
starting with the precooled molecule localized near the center of the trap, sideband 
cooling can then efficiently cool the molecule to the ground state. For the design of the 
optimal cooling strategies we should emphasize that the present setup does not only 
allow us to change A(t) and Q(t) in time, but also offers different possibilities to control 
7(i), A c (t) and n(t) during the cooling process. For example, using a tunable cavity 
[50| [5TI [52] the cavity frequency iv c (t) and therefore also the effective decay rate j(t) 
can be controlled in a time dependent way. In addition, when the cavity is coupled 
off-resonantly to another dissipative solid state device, e.g., a Cooper Pair Box, photon 
states will hybridize with the lossy two-level system and a substantial modulation of the 
effective photon loss rate can be achieved by changing the detuning between the two 
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Appendix A. Molecular master equation in the Lamb-Dicke limit 



In this section we present a few details of the Lamb-Dicke expansion (128]) of the molecular 
master equation ( |24l) which was derived in section [3j The first- and second-order terms 
in the Lamb-Dicke parameters i;,fj c <l read, 

Lip m = -I — 



a x {a + a)),p 



+ r] c g 2 (N + 1)\ (S_(z/)a + S_(-z/)a f )p m , a + + £_p m ,cr + (a + a 



+ H.c. 



+ Vc9 2 N^ (E+(V)a + S+(-^)a t )p m ,cr_ + £+p m , a_(a + a)) + #.c.j, (A.l 



and 



£ 2 p m = vl9 2 ( N + !n (£-0)a + S_(-z/)a f )p m ,cr + (a + a f ) 



H.c. 



+ Vc9 2 N 



+ H.c. )■ . 



(S + (z/)a + S + (-i/)a t )p m , a_(a + a 1 ") 
The operators S±(z^) are defined by, 

roo 

E ± (» = / rfre i(,/±(Ac - A))r e- KT a ± (-r), E± = E±(0) 



(A.2) 



(A.3) 



and depend on the Heisenberg operators c±(t) = e lHlt a±e~ lHlt with H\ = — (A/2)a z + 
(Q/2)a x the bare TLS Hamiltonian. For the evaluation of the integral ( 1A.3I) we follow 
reference [26] and write the exponential as 



' T ) 1 + isin (t 



A 



A A 



A = VA 2 + ft 2 . 



Using the algebra of Pauli matrices we arrive at the general form 
with the complex coefficients Bi{uj) given by 



(A.4) 



(A.5) 



4A 5 



2ft 2 



k — i(u — A r 



E- 



(A + sA) 5 



K-i(uj- (A g + sA)) 
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BoM 



B:Au) 



with A g = A c - A. 



1 



4A 2 
1 

AA 2 



2Q 2 



k — i(u — A r 



-2QA 



k — i(u — Ag) 



E 

s=± 

E 



« - z(u; - (A g + sA)) 

fi(A + sA) 
k - i(uj - (A g + sA)) 



(A.6) 



Appendix B. Bloch equations 

For later convenience, we write the Bloch equations (I3T1) corresponding to the effective 
TLS dynamics ( 1301) in the form 

(B.l) 



(a) = A(a) - T 



with a 



(a x , a y , a z ) T and 
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-(j N - lx )/2 


A + 


~5-6 x 


o ^ 
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r x 
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A = 




-(A + 5 + 5 y ) 


~(lN 


+ 7,)/2 


-fi 


f = 










V 


£l x 


n 


+ fly 


-In ) 




V 


7 


/ 



.(B.2) 



The parameters A and f2 describe the evolution of a driven free TLS. The interaction 
of the TLS with the cavity field introduces energy shifts 5, 5 X and S y , 



5 



-g 2 (2N + l)lm{B 1 (0)}, S x = 5 y = g 2 (2N + l)lm{£ 2 (0)}, (B.3) 



with the complex frequency dependent coefficients Bi(uS) presented in Appendix A 



Further, the effective decay via the cavity is described by the rates 7at = (2iV+l)7, J x , ly, 



7 = 2( ? 2 Re{ J B 1 (0)}, 



7a 



7, = 2< ? 2 (2A + l)Re{E 2 (0)}, 



and T x ,T y , 



T x = -2^Re{B 3 (0)}, T y = 2g 2 lm{B 3 (0)} . 

Finally there are additional effective Rabi frequencies fl x ,fl y , 

fl x = 2g 2 (2N + l)Re{5 3 (0)}, = 2<? 2 (2iV + l)Im{B 3 (0)}. 

From the definition of the Bi (1A.6I) we find that the non-standard terms j x ,y, 8x,y, flx,y 
and T Xt y vanish in the limit 11 C k or else for fl ~ k when the additional condition 
|A|, |A C |, I A c - A| > fl is fulfilled. 



(B.4) 
(B.5) 
(B.6) 



Appendix C. Evaluation of Sq(u) 

In this section we calculate the excitation spectrum defined in equation ([3 

v 2 fi 2 



Sq(uj) 



poo 

Re / drTrjKe^V^je^ 
Jo 
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This expression can be evaluated using the quantum regression theorem [49], which 
relates the two-point correlation function to the dynamics of the TLS represented by 
the Bloch equation (IB. If) . After performing the integral we obtain 

SnH = ^ Re ( . A , (C.l) 

v ; 2 \iu Det(iul + A) J ' v ; 

where the function in the numerator is given by, 

h{iu) = (1,0, 0)Ad(z^l + A) (iu(aa x ) + f<o- x > ) . (C.2) 

Here we also used the relation [icul + A}^ 1 = Ad(iul + A) /Det(iul + A), where Ad(x) 
is the adjugate matrix of x. By writing the determinant as 

Det(iwl + A) = (iu + e )(iu + e + )(iu + e_), (C.3) 

with €i the eigenvalues of A, we clearly see the three-resonance structure familiar from 
the excitation spectrum of a driven TLS. 

Simple analytic expressions for Sq(u) can be obtained in certain limits. In the weak 
driving regime Q <C 7/v, we can set Q = in the evaluation of h(iu) and and obtain the 
expression given in equation (|42j) . For general Q but in the limit 7jv A = VA 2 + Q 2 , 
the spectrum consists of three well separated peaks, with the center frequencies and 
the widths given by the imaginary and real parts of the 6j. For k such that the Bloch 
equation (IB.lj) acquires a standard form (see Appendix B ) we find 

e± ~ ± z'A - 7/2, 7 = 7tv (2 + sin 2 if) /2, 

e ~ iO - 70, 7o = Tat (l + cos 2 <p) /2, 

with simp = Q/A. Using a partial fraction decomposition of equation ( IC.lj) and keeping 
only the lowest relevant orders of 7 we end up with the simple structure for Sn(u>) given 
in equation (144j) . Since cooling and heating rates are determined by the peaks at non- 
zero frequencies, uj = ±A, we can approximately write the spectrum Sq(u) as 

r, 1 , an V 2 ^ 2 a ±l 2 fo (lilcosy?!) 2 \ AS 

Sn{ . K ±A) = V (. T aAw4 - 0± " C ° SV V' L + 2(2iV + l)(l + lv) J '< a4 » 
In section H] we use this expression to determine cooling and heating rates in the resolved 
sideband limit for arbitrary values of Q. 



Appendix D. Evaluate S g (uj) 

Let us now evaluate the V g spectrum (1381) , using the expression (1A.5I) for the generalized 
TLS operators E±(±i/). We write the resulting spectrum on the form, 

S g (u) = 2 V 2 c g 2 Re{(A+l) [B 1 (o;)p ee - B 3 (uj)p ge } + N [B x {-u) Pgg + B 3 (-iu)p eg } }, (D.l) 

where the complex functions Bi are defined in (1A.6I) . and the TLS steady-state 
expectation values pij = (i\Pi\j) are derived from the Bloch equations (j3l|) under the 
assumption Q k or Q, k |A|, |A c |. For simplicity we also assume 7/v A, which 
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is compatible with the constraints for CASC. For transparency we separate the thermal 
and the driven population, 



Pee/gg — Pee/gg ^ Pe 
N 



Pi 



N+l 



n 2 



2N + V P °" 2N+V Pee 2(2iV + l)(2A 2 + fi 2 )'^ D ' 2) 
Regarding the off-diagonal elements p e g/ge, it can be shown that only the real part 
contributes to the spectrum in the limit 7 at <C A, 



2A 

Re{p e9 } = Re{p ge } = — x p ee . 
Without further approximation, we obtain a spectrum with 6 resonances, 



(D.3) 
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vh 2 ^ 2 



2 A 2 
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kN 



k 2 + (lu 
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2A 2 

Vcf 

2A 2 



E K A +* A ) 
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E [( A + sA W - &p~] 



K + {UJ 
tt 2 p ee ] 



k(N - 



s=± 



k' + (u- (A g + sA)) 2 
kN 

K 2 + (u J + (A g + sA)f , 



(D.4) 



with A g = A c — A. The two limits of interest are the resolved sideband limit k < (/ 
and the Doppler limit 
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